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QQ ■ Abstract 

^P . Dark matter caustics are small scale, high density structures believed to exist in galaxies like 

ours. If the dark matter consists of Weakly Interacting Massive Particles, these caustics may be 

detected by means of the gamma rays produced by dark matter particle annihilation. We discuss 

C3 ' particle annihilation in outer and inner caustics and provide sky maps of the expected gamma ray 

distribution. 
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I. INTRODUCTION 

A consequence of the presence of cold dark matter in galactic lialos is the formation of 
caustics. Caustics are regions of infinite density in the limit where the dark matter parti- 
cles have zero velocity dispersion. The leading dark matter candidates, namely axions and 
Weakly Interacting Massive Particles(WIMPs), have very small primordial velocity disper- 
sions (~ 0.03 cm/s for a 100 GeV WIMP, ~ 10~^ cm/s for a 10 yueV axion [l|) and may 
be expected to form caustics with large density enhancements. If the dark matter is made 
up of WIMPs, the particles can annihilate producing gamma rays. It may be possible to 
detect this gamma ray flux by means of future experiments such as GLAST J2(]. Caustics 
may possibly be revealed by this technique. 

Caustics form because cold dark matter particles occupy a thin, three dimensional hy- 
persurface in phase space [l|, y] . There are two kinds of caustics in galactic halos: outer and 
inner jl| . The outer caustics are topological spheres surrounding galaxies. The inner caustics 
have a more complicated geometry which depends on the details of the dark matter angular 
momentum distribution jj^]. When the angular momentum distribution is dominated by a 
net rotational term, the inner caustics have the appearance of rings [l|, |j]. Observational 
evidence for caustic rings was found by analyzing the rotation curves of external galaxies 
|5j and the rotation curve of the Milky Way [6j. A triangular feature seen in the Infrared 
Astronomical Satellite (IRAS) map of the Milky Way is also suggestive of the presence of 
a nearby inner caustic ring [61]. Recently [3], possible evidence was found for a caustic ring 
in the galaxy cluster C10024+1654. In this article, we discuss dark matter annihilation in 
caustic spheres and caustic rings and present simulated gamma ray sky maps. 

The existence of dark matter caustics is sometimes disputed on the basis of large scale 
cosmological simulations. However, current simulations use particles with mass of order 
~ IO^Mq. Such massive particles have spurious collisions which destroy fine details in 
the dark matter phase space hypersurface. Thus the absence of caustics in most large scale 
simulations is probably due to the finite resolution of these simulations and the large particle 
masses involved. Note that discrete flows and caustics are observed in iV-body simulations 
when special care is taken. The simulations of Stiff and Widrow [8|], which increase the 
resolution in the relevant regions of phase space, do show discrete flows. More recently, 
Vogelsberger et al. [9| employed a new technique to study small scale structures such as 
caustics. Further progress in numerical techniques may enable the study of caustic properties 
in galactic halos using A^-body simulations. 

In cold dark matter cosmology, galactic halos form through a process of "hierarchical 
clustering". Small halos form first. Larger halos, such as that of the Milky Way, are the 
outcome of past mergers, the accretion of smaller halos of various sizes, and the accretion 
of unclustered dark matter. This general scenario must surely be correct without, however, 
telling us anything precise. The critical issue for the existence of a caustic is the velocity 
dispersion of the flow that would form it. To be specific, consider the flow of dark matter 
that the Milky Way halo is accreting today. If the dark matter in this flow has not clustered 
on any small scales yet, its velocity dispersion has the primordial value mentioned earlier, 
of order 0.03 cm/s for a 100 GeV WIMP. If it has clustered already, its effective velocity 
dispersion is the velocity dispersion of the structure it has become part of. The largest 
structures for which we have evidence of recent accretion are small sattelites, such as the 
Sagittarius dwarf, with velocity dispersion of order 10 km/s. So one may conservatively 
assume that the flow of WIMP dark matter accreting onto the Galaxy today has variable 



velocity dispersion between 0.03 cm/s and 10 km/s. 

Dark matter caustics are generically surfaces on one side of which the density diverges, 
in the hmit of zero velocity dispersion, as 

d{x) = -j= (1) 



where x is the distance to the surface and A is a constant called the "fold coefficient" . For 
finite velocity dispersion Sv, the dark matter caustics in a galactic halo are smoothed over 
a distance 

6x ~ —R (2) 

V 

where v is the fiow velocity (of order 300 km/s) and R is the overall size of the halo (of order 
300 kpc). The velocity dispersion cuts off the divergence of the density at the caustic, i.e. 
the -^ behaviour stops for x < 6x. For the above range of velocity dispersion, 6x ranges 
from 10^^ cm to 10 kpc. Since outer caustics have size of order hundreds of kpc, they are 
not erased by the expected velocity dispersion. Inner caustics which have size of order 10 
kpc are smoothed out only when the fiow is entirely composed of sattelites with velocity 
dispersion of order 10 km/s. 

As was mentioned, evidence for caustic rings of dark matter in the Milky Way was found 
in the existence and distribution of sharp rises in the Galactic rotation curve and the presence 
of a triangular feature in the IRAS map of the Galactic plane J6| . This evidence implies an 
upper limit on the velocity dispersion of order 50 m/s, or an upper limit on the smoothing 
scale 6x of order 10~^ i? ~ 30 pc. 

The WIMP annihilation rate per unit surface of caustic is proportional to the integral 

rfa;rf(a;)2^ / fL = AHn i^] . (3) 

JSx X \0XJ 

This equation shows that caustics enhance the annihilation signal only by a logarithmic 



factor [10| which varies from ~ 3 to 20, depending on the velocity dispersion. 

In this article, we limit ourselves to the simple cases of a spherically symmetric outer 
caustic and an axially symmetric inner caustic ring. Our symmetry assumptions are made 
solely for the purpose of computational ease. In previous work j^, we verified that caustics 
are stable under perturbations, in particular perturbations that break any assumed sym- 
metry. Moreover, the symmetry assumptions are consistent with the observational evidence 
mentioned earlier. 

Previous work on the possibility of detecting dark matter caustics by indirect means 
include [lO|, lul, 112, llS] . Caustics and cold fiows are also relevant to WIMP and axion direct 



detection experiments [li 

The number of photons received by a detector per unit area, per unit time, per unit solid 



angle and per unit energy due to WIMP annihilation is given by the expression! 111. |15 



{E^J,b) = -—(E^) X ^—^. (4) 
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The term dS/dE^ depends only on the particle physics: 



dS , . <crav> ^ dN^,^ 



dE^' " 2m\ Y dE^ 



< aaV > is the WIMP annihilation cross section times the relative velocity averaged over 
the momentum distribution of the WIMPs, m^ is the WIMP mass, dNj^i/dE^ is the number 
of photons produced per annihilation per unit energy E^ in the annihilation channel i, and 
bi is the branching fraction of channel i. 

The term E{l,b) is called the emission measure and depends solely on the dark matter 
distribution. 

E{l,b) = I dxd\x) (6) 

Jlos 

where x{l, b) is the distance measured along the line of sight(/os) in the direction of the 
galactic co-ordinates (/, b), and d is the dark matter density. Since all detectors have a finite 
resolution, the emission measure is often expressed in terms of < ii^ >, which represents an 
average over the resolution of the detector. In this article, we limit ourselves to computing 
E{l,b). 

As an example, let us compute the emission measure due to a smooth isotropic halo with 
the density profile 

d(r) = d{r,) i±^. (7) 

Consider an observer at a distance tq from the center of the halo, and a line of sight making 
an angle ip with the line joining the observer with the center (we choose ■i/' = to be the 
direction of the halo center), a is the core radius and keeps the density finite. We choose 
d = 0.3 GeV/cm^ when r = tq = 8.5 kpc. The emission measure E is plotted as a function 
of angle ip in Fig. [H for the cases a = 1, 2 and 4 kpc. 

II. OUTER CAUSTICS 

Outer caustics are topological spheres. They are fold catastrophes with the universal 
density profile 

d{r) = ^^e{R-r) (8) 

where d{r) is the dark matter density, r is the distance from the center of the sphere, 
measured along the radial direction, R is the radius of the sphere, A is a constant 
for the particular caustic, and B is the unit step function. For the self-similar infall 



model [l6|, [iTI, [iSj of the Milky Way halo, the radii of the first n = 1, 2, 3, ■ ■ ■ outer 
caustics Rn — {440, 260, 190, 150, 120, • ■ •} kpc [l8|] and the corresponding values of An 
~ {8, 9, 10, 10, 11, ■ ■ ■} X lO-^Mg/pc^/^ [i3. We will assume a thickness 6R for the caustic 
sphere, which sets a density cutoff equal to A/\''JR. 

A. Observer outside the caustic sphere 

Consider a spherical outer caustic of radius R and thickness 6R. Let an observer be 
located at a distance tq from the center of the sphere {tq > R). Consider a line of sight 
passing through the sphere making an angle tp with the line joining the observer with the 
center of the sphere (■?/' = is the direction of the center). The emission measure E is plotted 
in Fig. [2] for vq/R = 2, for the two cases 6R/R = 10~^ and 6R/R = 10~^. The maximum 



value of i? ~ 1.87A^JR/5R occurs when the line of sight is tangent to the sphere. Fig. [3] 
shows a simulated sky map. 



B. Observer inside the caustic sphere 

Let us now consider the case where the observer is within the caustic sphere. We once 
again assume a sphere of radius R and thickness SR. The observer is located at a distance 
To from the center (ro < R). As before, ip is the angle made by the line of sight with the line 
joining the observer with the center {ip = gives the direction of the center). Fig. H] shows 
the emission measure plotted as a function of angle ip for the three cases r^/R = 0.1,0.3 
and 0.6. Figs. [5] (a),(b),(c) show simulated skymaps. 

In both cases, the emission measure for the first few outer caustics of a halo like the 
Milky Way is quite small compared to that obtained from the smooth halo of Eq. [71 for 
6R/R > 10"^''. However, E would be larger for caustics of smaller radii and could possibly 
make a significant contribution to the total flux [13l |. 

III. INNER CAUSTICS 

The inner caustics are made up of sections of the higher order catastrophes. As men- 
tioned in the Introduction, the inner caustics resemble rings when the dark matter angular 
momentum distribution is dominated by a net rotational term. Here, we assume that this 
is the case. A caustic ring is a closed tube whose cross section has three cusps. Consider 
a caustic ring of radius a, horizontal extent p and vertical extent q. See [it] for details. 
We choose cylindrical co-ordinates {z,p = y/x^~+^). Let us define the two dimensionless 
variables X = {p — a)/p and Z = z/q. For simplicity, we assume axial symmetry about the 
Z axis and reflection symmetry about the Z = plane. With these assumptions, we can 
obtain an analytic expression for the density d close to a caustic ring, in the limit p,q '^ a. 
For the special case of Z = 0, the dark matter density is given by [l2| 

when X < 



rf(X,0) = ^| g(l + ;^) whenO<X<l (9) 

when X > 1. 



47rG 1-1^ V v^ 
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/ is a constant for the particular caustic ring. For the self-similar infall model of the Milky 
Way halo 0, [13, fn ~ {0.11,0.046,0.029,0.021,0.017,...} for the first few n = 1,2,3, 



dark matter flows. For the more general case of non-zero Z, the density (i(Z, X^) takes the 
form [l2| : 

where Tj are the real roots of the quartic equation 

97 

T^-2r3 + (i-x)r2-— ^2 = 0. (11) 

64 

Note that the assumption of reflection symmetry about the Z = plane implies that the 
dark matter density is a function of Z^. See Appendix |A] for a derivation of Eqs. [9]-[TT] 

A. Observer outside the caustic ring tube 

Consider an observer located at a distance tq from the center of a circular caustic ring 
with p = q and Tq > a + p. Fig. [6] shows simulated skymaps of the emission measure. The 



caustic ring is aligned with the galactic plane in (a) and inclined to the plane in (c). The 
emission measure is largest when the line of sight is tangent to the ring. The tricusp shape 
of the cross section is clearly visible in (b). 

B. Observer inside the caustic ring tube 

Let us now consider the case where the observer is located within the tricusp shaped tube 
of a circular caustic ring. The observer is in the plane of the ring at a distance ro from the 
center with a < tq < a + p. We also set p = q. Fig. [7] shows a simulated skymap of the 
emission measure. The bright ring and the line at 6 = occur when the lines of sight pass 
close to the cusps. 

We see that for the inner caustics, the emission measure can be much larger than that 
obtained from a smooth halo (except near the center of the halo). For a comparison of the 
annihilation flux with the measured gamma ray background, we refer the reader to |12l] . 

IV. CONCLUSION 

We have discussed dark matter annihilation in outer and inner caustics and provided 
simulated skymaps of the emission measure. In Section [TTl we considered spherical outer 
caustics. When the observer is located outside the sphere, the emission measure E is largest 
when the line of sight is tangent to the sphere, falling off rapidly to zero thereafter (Figs 
|2|3|) . When the observer is within the sphere, the emission measure only changes gradually 
with the angle ijj (Figs. Ii|5|) . 

In Section [nil we considered circular, symmetric tricusp ring caustics. When the observer 
is outside the caustic ring, the emission measure is largest when the line of sight is tangent 
to the ring. Looking tangentially, one may hope to identify the tricusp cross section of the 
ring[Fig. [6]^b)]. When the observer is within the tricusp, the emission measure is largest 
when the lines of sight pass close to the cusps, resulting in the bright ring seen in Fig. [71 We 
find that for inner caustics, E can be significantly larger than that produced by a smooth 
halo, especially in directions away from the halo center. 
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APPENDIX A: DARK MATTER DENSITY NEAR AN AXIALLY SYMMETRIC 
CAUSTIC RING. 

Consider a single cold flow of dark matter particles falling in and out of the halo, forming 
a caustic. Let us consider cylindrical co-ordinates (p, z) where p = \Jx^ + y^. For simplicity, 
we assume axial symmetry about the z axis and reflection symmetry about the z = plane. 



Let p = a be the caustic ring radius, defined as the point of closest approach of the particles 
in the z = plane (the particles with the most angular momentum are in this plane). 
Consider a reference sphere such that each particle in the flow passes through the sphere 
once. We can then assign to each particle in the flow, a two parameter label {a, r) which 
serves to identify the particle. We choose a = 7r/2 — 6 where 6 is the polar angle of the 
particle at the time it crossed the reference sphere, r = is the time when the particles 
just above the z = plane cross this plane. See [l[ for details. We may expand p and z in 
a Taylor series [l[, 



p{a,T) 



-u (r - To) - -sa 
bar 



(Al) 



where b,u,s,TQ and the ring radius a are constants for a given flow. In terms of these 
constants, we can express the horizontal (p) and vertical (g) extents of the tricusp: 
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Let us define the dimensionless variables X and Z 

p — a 



X 
Z 



P 



z 



(A3) 



Using Eq. IA2I and Eq. IA31 we may express Eq. lAll in terms of X and Z 

2 27 fbToa\ 



X 
Z 



(T-1 
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where T = t/tq. The density at points close to the caustic is given by [l| 

d{X, Z^ 



1 .^-^ dM cos a 

~p^ dVldt ID2I 



(A4) 



(A5) 



where \D2\ is the absolute value of the two dimensional Jacobian determinant 

\d{x,y)ld{a,T)\ 



\Dn 



2bp 1-X + 2T^ - 3T 



(A6) 



and the sum is over the individual dark matter flows that exist at that point. We use the 
self-similar infall halo model with angular momentum ITJ] to estimate the mass infall rate 



dM 
dVtdt 



fv 



"rot 

AttG' 



(A7) 



/ is a parameter that characterizes the density of the flow , v is the speed of the particles 
forming the caustic, and Vrot is the rotation velocity of the halo. 



1. Case 1: Z = 

Let us first consider the simple case of points in tlie Z = plane. Wlien Z = 0, tliere are 
two flows at each point when X > 1, four flows at each point for < X < 1 and two when 
X<0. 



a. X >1 



When X > 1, we have two possible solutions: T = 1 ± yX corresponding to the two 
flows that exist at each point. Using Eq lASj IA6I and |A7|, and approximating cos a ~ 1 and 
a ^ p we find the dark matter density at points close to the caustic 



d{X,0;X> 1] 



AtiG b ap X -1 



(A8) 



b. < X < 1 



In this range, we have four possible solutions: T = 0, 0, 1 ± yX. Summing over the four 
solutions, we find 



d(X,0;0<X< 1) 



A-nG b apl-X 



' 'i, ' 



X. 



(A9) 



c. X <0 



For X < 0, we once again have two solutions T = 0, 0. The dark matter density is: 

fVrot V 1 1 



d(X,0;X<0) 



47rG' b apl-X' 



(AlO) 



2. Case 2: Z ^ {) 



We now consider the general case of computing the dark matter density at points outside 
the Z = plane. For Z ^ 0, neither a nor T can be zero, so we may express a in terms of 
T in Eq. IA4[ to obtain a quartic equation for T: 



y4 _ 22-3 ^ Q _ ^^,^2 _ ^^ 
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0. 



The dark matter density is given by 

d{X, Z^) 



A-rG b 



4"2 1 



1 



X + 2Ti-?,Ti\ 



(All) 



(A12) 



where Tj are the real roots of Eq. lAlli The number of real roots is given by the sign of the 
discriminant 
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(A13) 



There are four real roots ii f3 > and two ii (3 < 0. 

To solve Eq. lAllt let us make the substitution Ti = T — 1/2, thereby eliminating the 
cubic term in Eq. lAllI 



Ti^ - ( X + - ) Ti^ - XTi - r = 



(AM) 



where T = -^Z'^ + -^ — -^. We may solve Eq. IA14I by expressing it as the product of two 



(A15) 



quadratics: 












(Ti^ + CTi + d) (Ti^ - CTi - M 


C^ solves the cubic 
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C^ - X2 = 



Let us define the two variables M and A^ 

M 



9^2 



3 / 16 

•x-^y_9z^/ 1 

3 / 16 V 2 



C^ has one real root 
2 



C^ = tix + -\ + (N + VN^-M^y^^ +[n~ VN^ - M3) 
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(A16) 
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if A^^ > M'^ and three real roots 



C'^-lX + -]+ 2VM cos 



2mT 1 
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(A20) 



for n = 0, 1, 2 if A^^ < M^ . D may be obtained using any of the values of C. The desired 
solution of the quartic Eq. lAllI is then 



T 



¥^±V(¥) -« 



^±i/(f +i 



(A21) 



Once the real roots Tj are known, Eq. IA12l and v ^ b may be used to obtain the dark matter 
density. 
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FIG. 1: Emission measure due to a smooth halo with the density profile of Eq. [7|with a = 1 kpc, 
2 kpc and 4 kpc. vq = 8.5 kpc and do = 0.3 GeV/cm^. 
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FIG. 2: Emission measure from a spherical outer caustic(observer outside the sphere) for SR/R 
IQ-^ and SR/R = W'^. ro/R = 2. The constant A was chosen to be 0.97 x lO-^^^^k^ 
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FIG. 3: Simulated skymap of the emission measure from a spherical outer caustic, with the observer 
being located outside the sphere. r^/R = 2,6R/R = W-^,A = 0.97 x IQ-^^^V^- The scale 
shows E I (GeV/cc)^ kpc. 
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FIG. 4: Emission measure from a spherical outer caustic (observer inside the sphere). Shown are 
three cases: r^jR = 0.1,0.3,0.6. The constant A was chosen to be 0.97 x 10^^— J^Vkpc and 
5R/R = 10-^ 
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(b) Tola = 0.3 
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FIG. 5: Simulated sky maps of the emission measure due to a spherical outer caustic, with the 
observer located within the sphere. Shown are the three cases: (a)ro/i? = 0.1, (6)ro/i? = 0.3 and 
(c)ro/i? = 0.6. A was chosen to be 0.97 x m~^Q^^/k^. 6R/R = 10"^. The scale shows E / 
(GeV/cc)2 kpc. 
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FIG. 6: Simulated sky maps of the emission measure from an inner caustic ring, with the observer 
located outside the ring. / was chosen to be 0.02 and p/q = 1. (a) shows a caustic ring aligned 
with the galactic plane, with a/ro = 15/17 and p/a = 1/15. The tricusp cross section of the ring 
is shown in detail in (b). (c) shows a caustic ring inclined to the galactic plane. The scale shows 
log,[-E/(GeV/cc)2 kpc]. 
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FIG. 7: Simulated sky map of the emission measure from an inner caustic ring, with the observer 
located within the tricusp tube. The ring is aligned with the galactic plane, p/q = 1 and / = 0.02. 
The emission measure is largest when the lines of sight pass close to the cusps, resulting in the 
bright ring seen in the figure. The scale shows logg[E/(GeV/cc)^ kpc]. 
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